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Abstract 

We consider a fluid queue where the input process consists of N identical sources 
that turn on and off at exponential waiting times. The server works at the constant 
rate c and an on source generates fluid at unit rate. This model was first formulated and 
analyzed by Anick, Mitra and Sondhi in We obtain an alternate representation of 
the joint steady state distribution of the buffer content and the number of on sources. 
This is given as a contour integral that we then analyze in the limit N — > oo. We 
give detailed asymptotic results for the joint distribution, as well as the associated 
marginal and conditional distributions. In particular, simple conditional limits laws 
are obtained. These shows how the buffer content behaves conditioned on the number 
of active sources and vice versa. Numerical comparisons show that our asymptotic 
results are very accurate even for N = 20. 



1 Introduction 



In traditional queueing theory one usually considers a system with one or more servers, at 
which customers arrive to receive some kind of service. As all servers may be busy, the 
customers may have to wait for service in a queue. Since there is uncertainty about the 
actual arrival times of the customers and/or their service requirements, queues are typically 
modeled stochastically, as random processes. 

In the last twenty years, models have appeared in which a continuous quantity, referred 
to as a fluid, takes the role of the number of customers or queue length. In these models, 
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fluid flows into a reservoir according to some stochastic process. The server may be viewed 
as a tap at the bottom of the reservoir, allowing fluid to flow out. The rate at which this 
happens is often constant, but may also be stochastic. Since the fluid reservoir takes the role 
of the traditional customer queue, it is often referred to as a fluid queue. A third term which 
is often encountered is fluid buffer, stressing the fact that the storage of fluid is temporary, 
to prevent loss of fluid at times when the input rate exceeds the output rate. Such models 
are used as approximations to discrete queuing models of manufacturing systems, high-speed 
data networks, transportation systems, in the theory of dams, etc. The distribution of the 
fluid level or buffer content, the average buffer content, the overflow probability (in case of a 
finite buffer) and the output process are the main quantities of interest in these fluid models. 

In the last decade, the literature on queueing theory has paid considerable attention to 
Markov-modulated fluid models jSO], |H]- I n these models, a fluid buffer is either filled or 
depleted, or both, at rates which are determined by the current state of a background Markov 
process, also called Markovian random environment. 



1.1 The Markov-modulated fluid model 

In this section we describe a general model of fluid entering and leaving a single buffer system 
[21]. Let X(t) denote the amount of fluid at time t in the buffer. Furthermore, let Z(t) be a 
continuous-time Markov process. Z(t) will be said to evolve "in the background". We will 
assume that Z(t) has a finite state space Kf with M = {0, 1, ... , iV}. 

The buffer content X(t) is regulated (or driven) by Z(t) in such a way that the net input 
rate into the buffer (i.e., the rate of change of its content) is d [Z(£)]. The function d(-) is 
called the drift function. When the buffer capacity is infinite, the dynamics of X(t) are given 
by 

dx r d[z(t)}, x(t)>o 

dt \ max{d[Z(t)],0}, X(t) = ' { } 

The condition at X(t) = ensures that the process X(t) does not become negative. When 
the buffer capacity is finite, say K, the dynamics are given by 




d[Z(t)}, 0<X(t)<K 
max {d [Z(t)\ , 0} , X(t) = 
min{d[Z(t)] , 0} , X(t) = K 



The condition at X(t) — K prevents the buffer content from exceeding K. 

One of the main reasons why Markov-modulated fluid models have attracted so much 
attention is that they are relevant for modelling telecommunications networks. In the litera- 
ture, more attention has been paid to models in which the buffer capacity is infinitely large, 
because this case is easier to analyze. For most practical situations in telecommunications 
the infinite buffer case is a good approximation for the finite buffer case, since overflow of the 
buffer is assumed to be extremely rare. For highly loss-sensitive traffic, allowed loss fractions 
in the order of 10~ 6 to 10~ 10 are typical [§]. 
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Unlike classical queueing models which assume renewal arrivals, Markov modulated fluid 
models can capture autocorrelations in arrival processes. The continuous nature of the fluid 
also makes them more tractable analytically. Many results have been obtained for a variety 
of fluid queueing systems 0, 0, [12!, [TH], [23], [23, |2Si, [HE], [HE], £3|, |Sj, HE]- 

In many applications, the process Z(t) evolves as a finite state birth-death process. Birth- 
death processes have wide applications in many practical problems [T^]. They can be re- 
garded as continuous time analogs of random walks. 

The paper which has become the main reference for birth-death fluid models was written 
by Anick, Mitra and Sondhi [3 . We shall refer to their model as the AMS model. The AMS 
model describes an infinitely large fluid buffer which is fed by N identical exponential on-off 
sources and emptied by an output channel with constant capacity c. An on source turns off 
at rate 1 and an off source turns on at rate A. Thus, the net input is regulated by a specific 
birth-death process Z(t) with state space Af, birth rate A& = A (N — k) , death rate //& = k 
and drift function d(k) = k — c. The rates are conditioned on Z(t) = k, < k < N. 

In a way, the AMS model was a generalization of the earlier model of Kosten rjU], |2Uj . 
where the limiting case N —>■ oo, A —>■ (with XN fixed) is considered. According to the 
authors in |3] , their model and its variants had also been proposed in other previous papers 

In [21], Kosten generalizes the AMS model by considering the same problem for a Multi 
Group Finite Source (MGFS) system, consisting of m groups of identical sources and a single 
shared buffer. In [22] Kosten considers a central processor, working at uniform speed, that 
receives information at uniform rates from a number of sources and a buffer which stores the 
information that cannot be handled directly. He shows that the stationary probability of 
overflow G(x) satisfies G{x) ~ Ce~ ax as x — > oo. For fairly complex aggregates of sources, a 
method is developed (called the decomposition method) to determine the exponent a. The 
coefficient C is determined by simulation. 

Daigle and Landford [TU] used the AMS model to study packet voice communication 
systems. Tucker jSU] also used a similar model, but the server capacity was defined as an 
integer number of information units and the model had finite buffer size. Results on the buffer 
content distribution, cell loss and delay were given, along with simulation comparisons. 

The feature common to the above models is that they all assume a continuous time 
distribution for the on-off intervals. On the other hand, Li [26 introduced a discrete time 
model with a finite number of on-off sources and geometric distributions for the on-off 
intervals. He assumed that in one time unit only one on-off source can change state. The 
channel capacity was assumed to be an integer number of sources and the buffer size could 
be either zero (burst switching-clipping case) or infinity (packet switching case). 

In [31] Mitra shows how to deal with states i £ M for which d(i) = 0. He also generalizes 
the AMS model in two ways. Again, he considers a buffer which receives input from N 
i.i.d. exponential sources, called producers. However, now the output is also assumed to be 
Markov- modulated, by M i.i.d. exponential consumers. A second generalization is that he 
considers the case of finite buffer capacity. 

In the manufacturing literature Wijngaard [55] has analyzed the model treated in [31] 
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for the case of one producing and one consuming machine. Sevast'yanov has given an 
innovative approximation for a production line model of many stages, which incorporates the 
exact solution of a pair of stages. All of this work assumes finite capacity buffers. Another 
extension of the AMS model is presented in where the authors allow the possibility of 
two sources of input, one slow source of a fluid and another source that generates packets. 

Some birth-death fluid queues with finite Af have also been proposed as approximative 
models, when the modulating process in the original model is a multidimensional Markov 
chain, see [21] • in [E] the authors allow the modulating process to be an arbitrary birth-death 
process on a finite state space, but require the drift function d(i) to have a particular sign 
structure. An analysis without any restriction on the sign structure of the drift vector has 
been given in the survey paper [51J. 

Until recently, authors have mostly focused on stationary probabilities. The transient 
analysis of stochastic fluid models, i.e., the analysis of the fluid distribution in the buffer at 
an arbitrary time t, is a complex and computationally intensive task. There are only a few 
papers dedicated to the transient behavior of fluid queues with Markov input. Narayan and 
Kulkarni [33] derive explicit expressions for the Laplace transform of the joint distribution 
of the first time the buffer becomes empty and the state of the Markov process at that time. 

The Laplace transform has been often used to evaluate the transient behavior of fluid flow 
models. In [HE] Ren and Kobayashi studied the transient distribution of the buffer content 
for the AMS model. The same authors deal with the case of multiple types of inputs in |3Tj . 
These studies have been extended to the Markov modulated input rate model by Tanaka, 
et. al., in [39] . 

Another approach to transient behavior is due to Sericola 02] • This leads to a numerically 
stable recursive method for computing transient and first passage distributions, and applies 
to general rate and drift functions. 

The equilibrium probability that the buffer content exceeds x, Pr [X(oo) > x] , in the 
AMS model has been studied asymptotically, as iV — > oo, by Weiss [53] and Morrison [32] . 
The exact expression in [3] is convenient for calculating this probability for moderate values 
of N. However, when N is large it is useful to have asymptotic approximations. 

Weiss (see also |46j ) analyzed the large N asymptotics using the theory of large 
deviations. This yields an approximation of the form Pr [X(oo) > NB] w exp [— 1(B) N] , 
where 1(B) is characterized as a solution of a variational problem. Overflow asymptotics of 
this type for more general Markov-modulated fluid models are derived in [2H]. Duffield JT] 
considers time-dependent asymptotics for the overflow probability. 

Morrison [32] investigated the AMS model starting from the exact expression in [3], 
using the Euler-Maclaurin formula. He constructed different asymptotic approximations to 
Pr [X(oo) > x] for the three scales x = 0(1/ N), x = 0(1) and x = O(N). 

The heavy traffic analysis of the AMS model, where N — ► oo and = + O (N^ 1 ^ 2 ) , 
was studied by Knessl and Morrison [IB] . In this limit, computing the joint steady-state 
distribution of the buffer content and the number of on sources can be reduced to solving a 
partial differential equation, which is forward parabolic in a part of the domain and backward 
parabolic in the remaining portion. 
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In this paper we analyze the joint distribution of the buffer content and number of on 
sources in the AMS model asymptotically, as the number of sources iV — > oo. Our results 
extend those obtained in |32j . 

In Section 2 we describe the AMS model in more detail and in Section 3 we derive an 
integral representation of the solution, from which the spectral representation obtained by 
the authors in [3] can be easily derived. In Sections 4-9 we obtain asymptotic expansions from 
the exact solution by using the saddle point method. We have to consider eleven relevant 
regions of the two-dimensional state space. In Section 10 we summarize our results and 
numerically compare our approximations with the exact solution. In Section 11 we analyze 
the conditional and marginal distributions. 



2 The Anick-Mitra-Sondhi (AMS) model 

In this model a data-handling switch receives messages from N mutually independent sources, 
which independently and asynchronously alternate between on and off states. The on and 
off periods are exponentially distributed for each source. These two distributions, while not 
necessarily identical, are common to all sources. 

The unit of time is selected to be the average on period and the average off period is 
denoted by j. Thus, if we denote by On (Off) the random variables measuring the elapsed 
time a source is in the on (off) state, we have 

T T 

Pr [On <T}= J e- s ds, Pr [Off < T] = J \e- Xs ds 
o o 

and 

E [Pr (a source is on)l = . (2) 

A + 1 

An on source will generate packets of information at a uniform rate that we assume, without 
loss of generality, to be 1 packet per unit of time. While a source is off, it generates no 
packets. Thus, when k sources are on simultaneously, the instantaneous receiving rate at 
the switch is k. The switch stores or buffers the incoming information that is in excess of the 
maximum transmission rate c of an output channel (thus, c is also the ratio of the output 
channel capacity to an on source's transmission rate). We assume that c is not an integer. 

Let Z(t) be the number of on sources at time t, and X(t) be the buffer content. As long 
as the buffer is not empty, the instantaneous rate of change of the buffer content is Z(t) — c. 
Once the buffer is empty, it remains so as long as Z(t) < c. We assume that the buffer 
capacity is infinite, so we need the following stability condition 

<7<1, 7 = ^7 (3) 



A + l N 
to ensure the existence of a steady state distribution. Note that c < N implies that the 
buffer may be non-empty and c > means that the output rate exceeds the average 

input rate. 
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Here we consider only the equilibrium probabilities, defined by 

F k (x) = lim Pr [Z(t) = k, X(t) < x] , < k < N, x > 0. 

t— >oo 

The balance equations are [3] 

^ ~ C) ~dt = X i N -( k - !)] F k-i + + 1) F k+l - [\{N -k) + k]F k , 0<k<N, (4) 
with 

F k (x) = 0, k<£[0,N]. (5) 

Moreover, if the number of sources on at any time exceeds c, then the buffer content increases 
and the buffer cannot stay empty. Hence, 

F k (0) =0, [c\ + f < k < N. (6) 

Also, 

*<»> = (TT5F0**. °- k - N (7) 

since .Ffc(oo) is the probability that k out of A" sources are on simultaneously and from (j2J) 
this follows a binomial distribution with parameter 

In |3| Anick, Mitra and Sondhi derived the spectral representation of the solution to 
(@J)-(II|). They wrote the equation in matrix form 

DF'(x) = M F(x), x > (8) 

where 

F(x) =[F (x),F 1 (x) } ... } F N (x)] T } 
D is an (N + 1) x (N + 1) diagonal matrix with entries 

Dy = (j -c)5y, 0<i,j<iV 

and M is an (AT + 1) x (N + I) tridiagonal matrix with entries 



My 



j, j = i + 1, < i < iV - I 
-[A(iV-j)+j], i = i, 0<z<iV 
A(iV-j), J = * - 1, l<^<iV 
0, |i-j|>2 



Using (jHJ), they represented the solution F(x) as 



7V-LcJ-l 

F(x)=F(oo)+ a ^ lX ^ ( 9 ) 



8=0 



6 



where Zi is a negative eigenvalue and 0» is the associated right eigenvector of the matrix 
D ~ l M 

They showed that and (jSJ imply conditions on the derivatives of Fn(x) at x — 0, 

rj n TP 

^f(0)=0, 0<n<N-[c\-l. (10) 

From (jl(J|) they found the coefficients a, in Q with the help of Vandermonde determinants. 
Asymptotic and numerical results for the overflow probability of a predefined buffer backlog 
were also presented. 



3 Exact solution 



Here we shall solve the problem flI|-(JZJ) by a different method than that used in j3J. This 
will lead to an alternate form of the solution that shall prove useful for asymptotic analysis. 
Looking for a solution to (@J) of the form 

F k (x) = e 6x h k (9) 

we get the following difference equation for h k {9) 

A (N - k + 1) h k -i + {k + l) h k+1 - [A (AT — k) + k + (k — c)6] h k = 0, < k < N. (11) 

To solve we represent h k {6) as an integral 



i fimw) 

c 



(12) 



where C is an small loop around the origin in the complex plane and H(9,w) is an analytic 
function of w inside C. Using ()12|) in and integrating by parts gives an ODE for H(6, w) 



dH (w + je-l- 0)N 
dw w 2 + (A — 9 — l)w — A 



H 



whose solution is 



with 



V{9) 
Ri(0) - 



H(9,w) = [1 - R 1 (9)w] NV ^ [1 + R 2 (6)w] N[1 ~ vm 
(27-1)0- A - 1" 



1 

' 2 

1 

2A 



A(0) 
[A(0) + (A - 1 - 



A(0) = a/(0 + 1 - A) 2 + 4A 



i? 2 (0) = — [A(9) - (A - 1 - 



(13) 

(14) 

(15) 
(16) 



Using the binomial theorem we can also represent hk{9) as 



(17) 



which shows that hk(9) is a polynomial in 9 of degree N — k, 
h N (9) = l, h N ^(9) = 1 + {1 ~ l)d N, 

^- a(6,) " N + N 



h N - k (9) = Cl (6)N + ■■■ + c^iO)^- 1 + [l + { \ k] ^ )e] N\ \<k<N. 

The polynomials h^iO) are related to the Krawtchouk polynomials |17j. The functions Ck{9) 
are entire (polynomial) functions of 9, that may be computed from (jl7j) . 
Since H is analytic inside C, we have 

h N +i(9) = — <pH(9,w)dw = 0. 
2m J 
c 

Imposing the boundary condition (jSJ) yields hk{9) — for k < — 1. Thus, 

<9 fc # 



(0,O) = O, k>N + l 



v(e j ) = +- r , o<j<N 



and therefore H(9,w) must be a polynomial in w of degree < N. In view of (fHj) we must 
have 

J_ 

iV' 

which gives the equation for the eigenvalues 9j. Since from (JJJ) we need Ff.(x) — * -Ffc(oo) for 
.t — > oo, we consider only the negative eigenvalues. Solving for 9j in (ITHj) we get 

^ = -a(-0<O, < j < iV - LcJ - 1 (19) 

where 

p + 2(A - 1)(1 - x)x + (1 - 2x) y/p 2 + 4A(1 - x)x 

crix) = ; — ^ 20 

V ; 2( 7 -x)(l- 7 -x) V ; 

and 

p = 7 - A + A 7 (21) 
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with < p < 1 from (J3j). Combining the above results and using (|7J), we find the spectral 
representation of the solution 

1 /N\ N ~ lchl 
F ^ = JTTW{k) Xk+ ^ "i'"-" 1 '!-^ ( 22 ) 

where the constants aj are to be determined. To find aj we employ arguments different from 
those in [3]. 

Taking the Laplace transform 



oo 

= J e^ x F k (x)dx 



in (£Q) and using the boundary conditions © for k > \_c\ + 1 gives 

(k - c)dT k = X[N-{k- 1)] + (k + l) T k+1 - [\{N - k) + k\ T k , \c\ + 1 < k < N. 

(23) 

Equation (|2~3*j) has two independent solutions, hk($) and h* k ({}), with given by (JT2J) and 

hlW = ±jH*{#,w)-± I dw (24) 
2m y W K+L 



c 



where C is a small loop around the origin and 

H*{^w) = [1 - A# 2 (tf)H w(e) [1 + Ai^H^ 1 ^ 1 . 
Using the binomial theorem we get from (J24j) 

KW = ^Wj (" l-^)]) [_*(,)]' [fllWr . (2S ) 
From (|25|) we have 

hm = \ [-A 2 + (7 A + A — 7) 1? - 7 #] iV + ^-^£iV 2 , 



= c *(tf)jV + ■ ■ • + c^i^JV*" 1 + (A N k , 1 < k < N 

where the c* k {d) are polynomials that can be identified from (|23j) using (|T5|) and f|16jl . There- 
fore, the general solution to (J23J) is given by 



J^(tf) = A^hi-ff) + A 2 (i?) h* k (&), [c\+l<k<N 
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with Axi^ff), A2 (1?) still undetermined. We note that hk and h* k are entire functions of d. 
Since h N+ i ($) = but h* N+1 ($) 7^ 0, we see from (JSJ) that we need A 2 ($) = 0, so that 



Taking the Laplace transform in (|22|) yields 

iV-|cJ-l 



(26) 



{l + X) N \kJ 



3=0 



(27) 



from which we see that the only singularities of J^ki'd) are simple poles at d = and 
*& = 6j, < j < N — |_cj — 1. Since /ifc(?9) is entire, Ai(i9) must have poles at these points 
also. We thus rewrite (l26l) as 



'jv-Lcj-i 

n ^ 



J=0 



where -8(0) is an entire function. Using the inversion formula, we can represent F k (x) in 
the form 

' N-[c]-l 



Br+ 



n — 

11 0,-0 



h k (<d)d$ 



{21 



where Br + is a vertical contour on which Re(0) > 0. From (|27jl the residue of J-k("&) at 
= is equal to 

(29) 



and from (|14|) we know that 
which implies that 



(l + X) N \k 
H(0,w) 

h k (0) = 



A 



A 



N 



k-N 



Thus, from (jzojl . we get 



5(0) 



A + 



A 



AT JV-LcJ-1 



(30) 



r) n e. 



3=0 



We determine by examining the limit •& — > 00, which correspond to considering the 

boundary condition at x = 0. 
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From (JT5j) and (fTo^) we can easily obtain the asymptotic expressions 



V(ti) ~ 7 -0-, i? -> oo 

17 

# l( tf)~I + (A-l)i tf^oo 
A A 17 



Setting u> = ^ in (JT2J) and using (pTTj) yields 



/lfc(l?) ~ 17 



AT-fc 



27TZ 



1 + 



A 



7V-fc+l 



We also have 



N-[c]-l 



n 



(l- 7 )iV 
N-Jfe 

'_l)JV-LcJ 



3=0 3 



li 



1? — > oo. 



1? — > oo. 



Combining and (pI3j) we get 



^ fc (#)~.B(#). 



A 



AT-fc 



Since ((BJ) implies that 
we obtain 



lim \dTk 

i9— >oo 



,1? L 

N — k 

o, L c J + 1 < fc < 



1? — > oo. 



S(i?) = o (i9 fc " LcJ ) , -> oo, L C J + 1 < fc < 

Setting fc = [cj + 1 we see that B{d) is an entire function that is o (i?) as it 1 — ► oo. 
generalized Liouville theorem |17| . £?(i?) must be a constant 



S(l7) = 5(0) 



A 



A + l 



n °t- 

3=0 



Therefore we can write 
F k (x) = 



N 



X + 1J 2m 



1? 



Br+ 



JV-|cJ-l 



n 

3=0 



e, 



6< -i? 



Closing Br + in the half-plane Re(i?) < we recover the spectral representation 



3=0 



with 



A \ nN ~^~ 1 Q. 



A + 



i n 



6j — 6j 

o * J 



In conclusion, we have proved the following theorem. 
Theorem 1 The solution of is given by 

"n-[c]-i 



F k (x) 



A \ N 1 



A + 1 J 2m J -d 

Br+ 



n 



Oj-0 
J=0 J 



h k (d)dd 



where Br + zs a vertical contour on which Re($) > 0, 

1 rH(#,w) 



2m J 



w 



N-k + l dw > 



(36) 



C is a small loop around the origin in the complex w— plane 

H($,w) = [1 - R 1 ($)w] NV W [1 + R 2 ($)w] N[1 - VW] 



1 + 



(2 7 -l)^-A-l 



, A(#) = v / (^ + l-A) 2 + 4A, 



«iW = ^ [A(0) + (A - 1 - 0)] , = ^ [A(^) - (A — 1 - i?)] 

^ = _a (jv) <0 ' o^^^-LcJ-i, 



and 



o~(x) 



p + 2(A - 1)(1 - z)z + (1 - 2x)^p 2 + 4A(1 -x) 



2(7 — x)(l — 7 — 2;) 

We have thus re-derived, using different arguments, the spectral representation 



F,(x) 



{l + X) N \k 



N 



A 



A + l 



N N-[c\-l 

E 

3=0 



N-[c\-l 



n 



$j — Bj 



e^ x h k {B ) (37) 



obtained by the authors in |3 a . We also derived the integral representation (|3l)Jl that will 
prove more useful than ()37j) for asymptotic analysis. We note that the coefficients in the 
sum in f!37j) alternate in sign and this makes it very difficult to obtain asymptotic results for 
Fk(x) from the spectral representation. 

We study asymptotic properties of F k (x) as N — > 00 for various ranges of k and x. We 
shall use (|36|. approximate the integrand for N large and then evaluate the integral over 
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Br + asymptotically. We will introduce the scaled variables (y, z) = (j^, -^j and consider the 
limit N — > oo, with y and z fixed. In Section 0] we will derive asymptotic approximations to 
Fk(x) that are valid in most of the strip 3) = {(y, z) : y > 0, < z < 1} , with the exception 

o 

of a few boundary regions and one curve in D (the interior of D). These other ranges are 
considered in Sections El - El 



4 The main approximation 

We first evaluate the product in for iV — > oo. 
Lemma 2 Let 



jv-Lcj-i 



0; 



n 0j -*- 

3=0 1 



Then, for N — > oo and a fixed > 9q, we have 

P(0) ~ P(0; AT) = 



-0n 



0-0, 



■exp [iV/x($)] 



^(l-2 7 ) + A + l ln 



2A(0) 



M0) = ~M7(l-7)0 + p] 

[ 7 (1- 7)0 + p] [A -1-0 + A 



(3? 



(39) 



(A - 1 - 0) p + (A + 1)' 7(1 - 7) + (1 - A) 7 (1 - 7)0 + A(0)5 



and 



0o 



7(1-7) 



<0, 5 = (1 -7) 2 A + 7 2 > 0, 
p = 7A + 7 - A. 



(40) 



VFe note from that p(0) = 0. 



Proof. Since we assumed that c is not an integer, we introduce the fractional part of c 
defined by 

a = c - [cj , < a < 1 (41) 
which allows us to write iV — |_cj as (1 — 7) AT + a. Since 



In [P(0)] = In 



(l-7)JV+a-l f k \ 

n ° \n) 

fc=o 17 ^ a VjvJ 



(l-7)iV+a-l 
fc=0 



(7 



(I) 
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we can apply the Euler-Maclaurin formula j3] to get 



In [P(tf)] = - In 



a{0) 



+ a (0) 



(l-7)AT+a-l 



+ 



hi 



+ In 



^(1-7 + ^) 
tf + (7 (l- 7 +^) 

dfc + o(l). 



Changing variables from /r to i = -| we obtain 



(l-7)AT+o!-l 



In 



(7 



(*) 



1— 7— (1— a)e 



dk 



In 



t? + a (x) 



dx 



with 



1 

£= N- 



If we denote the indefinite integral in (|42jl by 

= / In 

then we have, as e — > + , 



t? + a (x) 



1— 7— (1— a)e 




/ ln 


a(x) 


_$ + a(x)_ 










dx ~ I: [(1 - 7 )1 - (1 - a)Ii [(1 - 7 )1 £ - li(0), 



where Ii [(1 — 7) ] means lim li{z). From (|20|) we get 

2 t 1-7 



r(0)= J v tr(l- T -e)~^ e-0" 

7(1-7) £ 



with 

Therefore, 



7 + A — 7A. 



I'(l- 7 )= lim ln 

x 1 1— 7 



a (x) 



$ + a(x) 



ln(l) = 0. 



To find li(z), we use integration by parts 



In 



a (x) 



+ 1 



dx = —z ln 



r) 



(TiZ 



d I x 



a \x\ 



\Q + a (x)} a (x 



-dx. 
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Since (fTHj) implies that V [—a(x)] = x, we can write (jlHj) as 

V [-a(x)} 



li{z) = —z In 















\& + cr (x)\ a (x 



-a' [x) dx 



or, changing the integrating variable from x to a, 



li{z) = —zln 



—z In 



a{z) 



V{-a) 



da 







a(z) 




■d 








f 1 




i 


[a(z) 









-zln 












a{z) 





(2 7 -l)a + A + l 
+ 1 - A) 2 + 4A 

a(z) 

— : — da 



($ + a)a 



da 



(49) 



+ 



■d f (2 7 -l)cr + A+l 



2 J y/(-a + 1 - A) 2 + 4A (0 + <r) 



The first integral in (|49j) is 



-dcr = — In 



cr 



(0 + a) cr d 

The second integral in ()49j) is quite complicated and after some calculation we find that 

(2 7 -l)cr + A+l 1 



a/(-(J + 1 - A) 2 + 4A + cr) cr 



d(7 



2$ 

fl(l-2 7 ) + A + l 
tfA(tf) 



1 , f [(A + l) 2 + (A-l)cr-(A + l)A(- 
m ' 



-a 



In 



4A(7 2 

, (A - 1 - g) cr + (A + l) 2 + (1 - A) g + A(fl) A(-cr) 



Hence, 





In 


r ] 


K 




j& + a{z)_ 







+ Qi Hz)] 



with 



n( s 1, /[(A + l) 2 + (A-l)a-(A + l)A(- 
Qi{a) = - In J 



a 



4 

+ 2A(tf) ln 



4Aa 2 

(A - 1 - 0) cr + (A + l) 2 + (1 - A) ■& + A(i?) A(-a) 

tf + CT 



(50) 



(51) 
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As z t 1 — 7 we have 

Qi [*(*)] ~ -~ ^ (A) + ^ (1 ~ 2 2 2|^ A + 1 ^ {2 [A - 1 - * + A(0)]} 
and taking (|37j) into account we get 

Ii [(1 " 7)1 = -~ ^ (A) + ^ {1 ~2A^) X + 1 ln {2 [A - 1 - ^ + A(0)]} . 
Using in (joTj) we have 



Q 1 W0)) = \^K 



2X 0(1 -2 7 ) + A + 1 



2A(#) 



x In 



, (A - 1 - g) p + (A + l) 2 7 (1 - 7) + (1 - A) 7 (1 - 7)1? + A(&)6 

7(1 - 7)1? + p 



From (03) and (|SUJ) we get 



li(0) = --ln 



P 



7(1 - 7)1? + p 



+ Qi b(o)] 



1, r v. q 1 1, ... tf(l-2 7 ) + A + l 
-ln[ 7 (l-7M + p]- i ln(A)+ 



x In 



, (A - 1 - g) p + (A + l) 2 7(1 - 7) + (1 - A) 7(1 - 7 )fl + A(g)j 

7(1 - 7)1? + p 



Combining these results in yields 



1— 7— (1— a)e 




/ h 









x In 



2 2 ^ A + 1 ln{2[A-l-^ + AW]} 

1, r M v. q , tf(l-2 7 ) + A + l 
--ln[ 7 (l-7)^ + p]- 2A ^ 

, (A - 1 - ■&) p + (A + l) 2 7(1 - 7) + (1 - A) 7(1 - j)# + A(#)5 

7(1- 7)1? + p 



1, , Na . tf(l-2 7 )+A + l 

- -In [7(1-7)1? + P\ + 



x In 



2 L,v ,y rj 2A(0) 

[ 7 (1~7)^ + P] [A - 1 - + AQ9)] " 

(A - 1 - 1?) p + (A + l) 2 7(1 - 7) + (1 - A) 7(1 - 7)1? + A(0)<f_ 
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From (J35j) we see that as iV — > oo 



1 
2 



In 



a(0) 



.i? + «t(0) 
Therefore, we conclude that 



In [P(i?)] = - In 



+ In 



^(1-7 + ^) 
tf + a(l- 7 + ^) 



In 



P 



_7(1 - 7)1? + p 



P 



1 , r m no 1 l^(l-2 7 )+A+l 

-^In[7(l-TM + P] + - 



x In 



_7(1 - 7)1? + p 

[7(1 - 7)1? + p] [A - 1 - + A 



(A - 1 - ■&) p+ (A + iy 7(1 - 7) + (1 - A) 7(1 - 7)0 + A W 



Exponentiating this yields the formula for P($; iV). 
We note from (JHHJ) that 



2 
5' 



and that P($; iV) has a branch point singularity at $ = # - 

We evaluate next in (|36|1 for N ^ 00 and = 0(N) with 1? fixed. We introduce 

the new variables 

x k 



y = —, z = —, y>0, < z<\ 
N N ~ ~ ~ 

and consider the asymptotic approximation to Fk(x) inside the domain 

£> = {(y,z) : y > 0, 0<2;<1}. 
Writing h^ft) in terms of z we have 



(52) 



(53) 



2m J 



w 



N{l-z)+l 



dw 



1 



w 



exp [A^?7(u7, 1?, z)] dw 



(54) 



with 



r/(w, 0, z) = V($) In [1 - Pi(i9)w] + [1 - V(0)] In [1 + P 2 (0)w] - (1 - z) ln(iu) 



(55) 



We first assume that 7 < z < 1 and compute the integral (J53)l as iV — ► 00 by the saddle 
point method [Hj. Thus, we locate the saddle points W(0, z) by solving the equation 



d_ 

dw 



rj(w, 1?, z) 



0. 



or 



zW 2 + [(7 - z) •& + z\ - z - A] W + (1 - z) A = 



(56) 
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and find two saddle points W + z) and W- (■#, z) denned by 

W ± (0, z) = ^ + ^~ A + ±- \\ - ^ ± y/W^j] (57) 



2 2z 

with 



2) = p 2 + [2 (A + 1) p + 20tf] - 7 ) + [(A + l) 2 + 2 (1 - A) ■& + tf 2 ] (z - 7 ) 2 . (58) 

We next study the motion of the saddle points W± (i?, z) , and also the branch points 
~R^W) anc ^ WTff) °^ r /( u '''^' z )' as ^ decreases from oo to 9 , with z > 7. First we observe 
from (|TH|) that for all real 0, Ri("d) and #2(0) > and thus 



1 <o< 1 



For large 1? we have 

iti z z — 7 17 

and hence for all 1? > 

-— < < < W + < -j-. (59) 
When 1? = the branch point coalesces with the saddle point W + 

-1- = 1, W + (0,z) = l, W-(p,z) = -{l-z) 
Ri(0) z 

and we get 

-RW) <0<w ' i0 - z)<w+i0 ' z) = Mf)- m 

For #0 < 1? < 0, the saddle point W+ moves to the right of the branch point "^^y 

< < < 4- < W+. (61) 
Ri R\ 

Finally, when d — 9q with 7 < z < ^- the branch point coalesces with the saddle point 
W- 

7 t 7 Vt = -(1-7), W + (0o,z)= i 7 ^^, WL(*o, *) = -(! -7) 
i2i(0 o ) 7 (1-7)* 7 

and we have 

-pTflT < < = ttTyTT < W+Cflo,*)- (62) 
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Note that when i? = where 



A(»-i)-»+2ya^y < y 

z — 7 (3 
the discriminant D($,z) vanishes and the saddle points coalesce. Thus, 

1 



(63) 



< < 



1 



Furthermore, 



since 



an 

+ A (227 — z — 7) 



(64) 



<9z 



2 — 7 



y/Xzll-z) (2- 7 ) 2 
0^(1-2) + A (227 - 2 - 7) 



Hence, D($, z) has a double zero at 1? = •&*. 

From (|59|) - (}62|) we conclude that for $ > 6q with 2 > 7 we must deform the contour C to 
the steepest descent contour through the saddle W- and therefore 



h k (0) 



exp [N7}(W-,ti,z)] 



d 2 rj 




dw 2 


w=W- 



■& > 6 , 7 < z < 1. 



(65) 



Note that the steepest descent directions at w = W_ are arg (w — W^) = ±| and we can 
shift C into another circle that goes through this saddle point. 
Using (|55|) and (|56j) we find that 



<9-u; 2 



A(ti,z)W± + B('&, z) 
(W ± f [(t? 7 - A) W± + A] 2 



(66) 



with 



2) = \ 2 {z - 1) (zA + z - A) + [(-2A 2 + A - 7A 2 - i)z 2 + A (4 7 A + 2A - 7) z - 3 7 A 2 ] 1? 
+ (2 - 7) [(27 A + A - 2 7 ) z - 37A] tf 2 - 7 (z - 7) 2 tf 3 



, z) = A(l - 2) { (zA + 2 - A) A + [2 7 A + (7 - A - 7A) z\ ■& + 7 (z - 7) tf 2 } . 
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Using and in we have 



A 



N 



A + ly y/2^N2iri J $W- V # - 0o 



1 1 



-0o 



Br" 1 



d 2 i] 



dw 2 



w=W- 



(67) 



x exp {N [ytl + + rj(W-, z)\}d&, z > 7. 



To compute the above integral as N — ► 00, we shall again use the saddle point method. We 
obtain Fk(x) ~ G\(y,z) with 



Gi(y,z) 



1 



A 



N 



1 



1 



2-nN \X + lJ QW4<3),z) V @ - 6 \j r/ ww (W_,Q, z) 
1 

■exp[JV*(y,W_,e,z)], 6 > 0, z>7 



(68) 



where 

w, 1?, z) = yd + //(i?) + 77 (it;, 1?, z) 
and 0(?/, z) is defined to be the solution of the equation 

3n 



(69) 



(70) 



where we have used 



dr] 
dw 



[W.(Q,z),Q,z] = 0. 



For each value of G, f)7Uj) defines implicitly a curve in the (y, z)-plane. 

Since (JBTj) has a pole at 1? = 0, (jSHj) is only valid if y and z are such that > 0, i.e., 
the saddle lies to the right of the pole. Then Br + can be deformed to the steepest descent 
contour through •§ = 0. We also note that when = 0, f!68j) is singular. We shall find an 
asymptotic approximation valid when the saddle is close to = in Section El From (J70*|) 
we see that = corresponds to the curve 



y = Y (z) 



z — 7 



P 



A + l (A + l)' 



In 



zX + z — A 
P 



7 < z < 1. 



(71) 



For ?/ > 5^0(2) the saddle point lies on the negative real axis, to the left of the pole at 
•d = 0. We re- write (JHEJ) as 



F k (x) - 



l + X) N \k 



N 



A 



A 



A + l 



iV 



2vu 



x) 



Br" 



JV-[cJ-l 



n 



j=o J 



h k {&)dti (72) 
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where #0 < Re($) < on Br~. Note that the residue at $ = was already calculated in 
(J2HJ) and corresponds to F fc (oo). Using (|55jl and (|6*5|) in (j72Jl and approximating the integral 
by the saddle point method, we obtain 



F k (x) - 



(l + X) N \k 



N 



A 



N 



-On 



1 



2nN \x + ij e w_(e, z) v 6 - 0o y WW-> e, z) 

-exp[JV*(y,W_,e,2)], O <©<O, 7<z<l 



(73) 



with ^(y, iu, 1?, z) and ©(?/, z) defined by ((OHj) and (ffHjl respectively. 
From (fTU|) we find that the value = 6q corresponds to the curve 



v , \ [7(1 -l)f P \ 1 2 -Sz 
y = Yi(z) = { — + In 



S 2 



7(1-7)/) 



7(1 -l)p 

7 2 — 5z 



-I}. 7 < z < j • (74) 



Our approximation to the integral in (136(1 assumed that $ > 6q. However, ()73)1 passes 
smoothly through the curve y = Y\(z) since from we have 



d 2 n 7 



2 \„,2 j_x2 (T 



;i- 7 ) z A7 2 + 5 



— z 



A 2(l _ 7 )2^2 (tf _0 Q ) 



tf-+0o. 



The discriminant Z?(0, z) vanishes when 



drj 



A (1 — z) 



,r,z 



7 

7 < ^ < — 



(75) 



where $*(z) was defined in (|U3*j) . But, as we showed in (j64|). the derivative of the discriminant 
also vanishes along the curve y = Y*(z) so that (|75|) can be smoothly continued through this 
curve by simply replacing W_ by W + . 
We summarize our results below. 

o 

Theorem 3 Defining y — j| and z = ^, /or (y, z) G D = {(y, z) : y > 0, < 2 < 1} we 
have 



1. 



F k {x)r,G 1 (y, z) 



A \ N 1 



-0n 



V M (y,W-,Q,z 



2nN \X + lJ W_( e > z) V © - O V ^(W-, 0, z) 
exp[iV*(y,W_,0,z)], (y, z) £ 9*i 



(76) 
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1-1 


0>O 


Yn(z) / 

0=0/ 0<O 






0.8- 






w. 




0.6 






W+ = w. 


Y*(z) 


z 










0.4- 






W+ 




0.2- 















0.2 0.4 0.6 

y 


0.8 


1 



Figure 1: A sketch of the regions 9ti, 9^2 and 9^3 



Xi = {(y,z):0<y<Y o (z), 7 < 2 < 1} , 



(77) 



^-A + l (A + l) 2 

where 

7](w, tf, *) = In [1 - i2i(#H + [1 - V(ti)] In [1 + R 2 {$)w] - (1 - z) \n(w) 

•d + 1 - A 1 



W±(t?,z) 



2 + 2z 



A - 71? ± a/^M 



£>(0, z) = p 2 + [2 (A + 1) p + 20tf] (z - 7) + [(A + l) 2 + 2 (1 - A) # + tf 2 ] (2 - 7 ) 2 . 

w, z) — y& + + 77 (w, 2) 
and 0(y, zs defined to be the solution of the equation 

y + fi'(Q) + ^[W-(e,z),e,z] = 
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^(l-2 7 ) + A + l 1 
+ 2M¥) M 



^) = --ln [ 7 (l- 7 )^ + p ] 

[ 7 (1 - jy& + p] [A - 1 - g + AQ9)] " 

L(A - 1 - 1?) p + (A + l) 2 7(1 - 7) + (1 - A) 7(1 - 7)1? + A(tf)5_ 



A fc ~Gi(3/,2), (y,z) G9* 2 



9l 2 = {(y, *) : Y {z) <y<Y*(z), 7 < * < 1} 



with 



Y*(z) = -im'(&*) 



dr] 



X(l-z) 



7 < z < 



r 



3. 



F k (x) 



X(z-l)-z + 2y/zX(l-z) 



z — 7 



7 < z < 



7 



(y, *) : 2/ > Y*(s), 7 < 2 < y ^ U {( ij. : ) : < //. <' : < " | 



where 



G 2 (y,z) 



A 



A 1 



-00 



2vrAr va + 1; 0+ w+ (e+, *) v e+ - e v ^(W-h e+, *) 
ex P [m(y,w+,e + , 



and Q + (y,z) is defined to be the solution of the equation 

y + tl '(e + ) + ^[w + (e + ,z),e\z]=o. 



(78) 
(79) 



(80) 
(81) 

(82) 



(83) 



The three regions %Kj are sketched in Figure ^ Note that as we pass from 9ti to 9K 2 there 
is a "phase transition" in the asymptotics. However, the solution is smooth as we pass from 
9*2 to 9t 3 . 

Theorem |3] ceases to be valid along the curve y = Y (z), 7 < z < 1, because the 
function Gi(y, z) is singular there. We shall find a corner layer solution near the point (0, 7) 
in Section a transition layer solution along y = Y (z), 7 < z < 1 in Section and a 
corner layer solution near the point (Yn(l), 1) in Section I5~2l 
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From (|57j) we have 



W-(#,z)~ 1 + ^ (l-z), ^>^o (84) 

w+ (0, z) '—, z^O, $<9 Q 



and from (J55J) we get 



[1 + (1-7)*1 



2 



^(wi,^) ~ 1 v ,y J , $>e (85) 

A z 1 — z) 



Hence 



r] ww (W +} $ } z) ~ — -2, 2^0, ^ < 6> . 

(A - 71?) 



G 1 (y,z) = 0[-^=Y z^l 
G a (y,z) = 0[±=Y z^O 



and our approximation develops singularities as z — > and z — > 1, which corresponds to 
fc ~ and ~ iV. We shall find a boundary layer solution near z = in Section and a 
boundary layer solution near z — 1 in Section |HJ In these cases we will need to re-examine 
the expansion of 

Finally, from (J70j) we see that 0— » oo as y — > if j < z < 1. We shall find a boundary 
layer approximation near y = with 2 > 7 in Section El To analyze the case of small y 
we will need an expansion of P(i?) in Lemma El valid for d large. Note that Theorem El 
does apply for y — > as long as < z < 7. In this range we can show that G + in (|83p 
remains finite. In particular, the boundary masses F k (0) for k = Nz with < z < 7 can be 
computed by setting y = in (J5U|) . We also note that in this range F k (0) ~ F k (oo) and the 
exponentially small (as N — > 00) difference can be estimated from (|8()J) . 

5 The corner layer at (0,7) (Region I) 



We next consider the vicinity of the corner point (y,z) = (0,7). We first evaluate the 
product P(tf) for iV -> 00 and ■& = 0(N). 



Lemma 4 Let $ = JVS, 5 > 0, 5 = 0(1) and 



tf-lcj-i « 

p(jvs) = = n 



i=o J 



04-NS' 
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Then P(S) ~ P(S; N) as N -> oo, w/iere 
1 



07d-7) V^ 1 -P{"^ d-7) m [(1-7)5^-^^7)} 



x r ( — + 1 — a) exp < — In 



S 



s 



7 



L0(1-7)7V 



+ 



20 - p - 1 
5 



i/ere Y (•) zs i/ie Gamma function, a = c — [c\ is the fractional part of c, <p = 7 + A 
and p = 7A + 7 — A. 



Proof. We write 



(l-7)7V+a-l 



r) ■ 



i=o J 



(l-7)JV+a-l 

n 

j=0 



-a 



a) 



-a (£) - SN 



(l-7)JV+a-l 



n 



(7 



(£) (1-7-^) 



|i [(l-7)tf-j]S + a(£)(l-7-£)' 
Changing index from j to n = (1 — 7)iV + a — 1 — j and defining 

a — 1 — n 



a n = a 1 - 7 + 



we have 



(l-7)iV+a-l 

= n 



iV 

V N ) 



n=0 



(l-7)AT+a-l 



n 



n=0 
(l-7)JV+a-l 



1 / n + 1 — a 



n (n+ l_ a)5+ (nto) 



+ 1 - a) 5 + (2^) a n 

(l-7)AT+a-l 

— a \ 

cr n X 

n + 1 — a)S 



n 



n=0 



+ (n + 1 - a)S 



Or. 



V 1 (S)P 2 (S)P 3 (S). 



We evaluate individually the three products in (jHHj) beginning with Vi(S). 
Setting e = -k and using the Euler-MacLaurin summation formula, we have 



In 



(l-7)jV+a-l 

] | e (n + 1 — a) a n 

n=0 



(l-7)JV+a-l 

2j In [e (n + 1 — a) a n ] 

n=0 



- In [(1 - 7) a (1 _ 7)iV + a -i] + - In [e (1 - a) <t ] 

(l-7)JV+a-l 



+ J ln{e(n + l-a)<r[l-7 + (a-l-n)e]}dn + o(l). 
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Changing variables to x = 1 — 7 + (a — 1 — n) e , we get 

(l~7)iV+a-l 



J In {e (n + 1 — a) a [1 — 7 + (a — 1 — n) e]} dn 





l-7+(a-l)e 



1 

e 



In [(1 — 7 — x) a (x)] dx. 



z 

hbz) = J hi[(l — j — x)a(x))dx, 



Defining 
we have 

l-7+(a-l)g 

J In [(1 - 7 - x) o {x)\ dx ~ l 2 (1 - 7) - (1 - «)I' 2 (1 - 7) £ ~ I 2 (0), 


From (|20j) we have 

, , , p + 2(A - 1)(1 - x)z + (1 - 2z) ^p 2 + 4A(1 - x)x 
(l-l-x)a(x) = ^—^ 

and evaluating this at x = 1 - 7 we get 

r 2 (l- 7 )=ln(0). 

After some calculations, we find that 

h{z) = (7-*)ln(2)-z-ln(l-z) 

+ I In [p 2 - 2A (2^7 - 7 - *) + ^(z)] 

+ -(7 - 1) In [p 2 + 2A (2z 7 - 7 - z + 1) + ^(z)] 



+ 2; In 



p + 2z (1 - z) (A - 1) + (1 - 2z)K 1 (z) 

7 — z 



+ Un[p 2 + 2\(l-z)+pK 1 (z)] 

where 

K x (z) = ^/p 2 + A\z(l-z). 

From (|89|) we obtain 

I 2 (0) = 2 7 ln(2) + ha + 1) MA + 1) + I In (5) + 7 ln( 7 ) 
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and 



I 2 (l - 7) = 2 7 ln(2) + 1( 7 + 1) ln(A + 1) + 1 In (5) + 7 - 1. 



Thus, 



I 2 (l- 7 )-I 2 (0) = 7 -l- 7 ln( 7 ) 



and we have 

l-7+(a-l)e 



If 1 

- / ln[(l - 7 - a;)cr (a;)] da: ~ [7 - 1 - 7 ln( 7 )] - + (a - 1) ln(0), e -> 0. 



From (HO]) and (JHZ|) we get 



- {In [(1 - 7) a ( i„ 7)A r +Q _i] + In [(1 - a) ea ] } 
= iln[(l- 7 ) (7 (0)] + iln{(l-a)ea[l-7 + (a-l)e]}~^ln 



7 



Therefore we conclude that 

~(l-7)7V+a-l 



hi 



n 

n=0 



n + 1 — a 
JV 



On 



[7-I-7 ln( 7 )] JV + {a - 1) ln(0) + - In 



7 



e -> 0. 



and then 



(l-7)AT+a-l 



n 

n=0 



1 / n + 1 — cc 

~lv 



(l-7)JV+a-l 



-(l-*/)N-a 



n 

n=0 



iV -> 00 



n + 1 — a 



JV 



On (90) 



— exp [-(1 - 7) ln(0) JV + (7 - 1)JV - 7 \n0y)N] , JV -> 00. 
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We next consider the second product in (|88j). Using properties of Gamma functions [1 
we have 



(l-f)N+a-l 

W) = n 



n=0 



+ ^(ti + I -a) 



(l-7)JV+a 



r(| + i- a) 
r[(i- T )Jv+f + i] 



and Stirling's formula [Mj gives 

1 



V2nN \S 



S 



5) r(^ + i-«)[(i- 7 )iv]^ 1 



(91) 



x exp <j (1 - 7) In ( I ) JV + (1 - 7)JV - (1 - 7) [In (1 - 7) JV] 



27 



Finally we determine the asymptotic approximation to V^S). The Euler-MacLaurin sum- 
mation formula yields 



(l-7)iV+a-l 



n=0 



<f)+(n+l-a)S 
(n + 1 - a) (S + ea n ) 



ln 



0+(l- 7 )S7V 



(1 - 7) (SN + (T(i_ 7 )Ar +a _i) 



In 



+ / >„ 





+ (n + 1 - a)S 



(n + 1 - a) (5 + so [1 - 7 + (a - 1 - n) e]) 



0+(l-a)5 

(1 - a) (S + £<7 ). 

c?n + o(l). 



(92) 



Changing variables in the integral in (}9~2*j) from n to x = 1 — 7 + (a — 1 — n) e, we have 

(l-7)AT+a-l 



hi 



+ (n + 1 - a)5 



n + 1 - a) (£ + eo" [1 - 7 + (a - 1 - n) e]) 

l-7+(a-l)e 



dn 



hi 



+ S(l - 7 - x) 



(1 - 7 - x) [etr(x) + 5] 


We write (|9*3*)) as the difference of two integrals 

l-7+(a-l)e 



dx. 



l-7+(a-l)e 



hi 



0£+ (1 -7-x)S 
(1 - 7 - x) [£(j(x) + S] 

l-7+(a-l)e 



G?X 



1 — 7 — x 



+ 5 



dx 



ln [£cr(x) + S] dx. 




The first integral in (jHljl is easily evaluated as £ — ► to give 



l-7+(a-l)e 



hi 



0, f (l- 7 )e5 



+ 5 



1 — 7 — x 

£ + (1 — a) ln 



dac = (1 - 7) In (S) 
1 — a 



Introducing the function 

z 

I 3 (z) = / In [ea{x) + S] dx, 



£ + 0(£ 



(93) 



(94) 
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in (|9~3jl and using (fT3j) we have 

z 

h(z) = z\n [ea(z) + S]- J 



ex 



ea{x) + S 



a'(x)dx 



a(z) 



zln [ea{z) + 5] 



ea + S 



da. 



After some work we find that 



/ 1 \ q r / >, 01 (27-l)5-(A + l)e 
I 3 (2;) = [z--)\n [ea(z) + S] + 



x ln< 2 



2Q 2 

A [-a(z)\ Q 2 + [(A — 1) e — 5] a(z) + (1 - A) 5 + (A + l) 2 e 



+ 7 



9 



ea(z) + 5 
\n{X-l+a(z)+A[-a(z))} 



with 



Q2{3) 

From (|95j) we have, as e — > 



[(A + l)e] 2 -2(A- l)eS + S 2 . 



I 3 [l- 7 +(a_l) e ] 



7 



In 



iA\, <p Ue\[(j)+(l-a)S} 

~ ~ 11 05 



and 



ls(0) 



+ (1 - a) In 



5 



7 I In I — I — ^— In 



X S 



S J s 

1 — a 
+ (1 -a)5 

"4A (1 - 7) 

7 



£+ (2 7 -l)(A-l) £ 



£ + 



5 (A-D-P 



5' 



We conclude that 

l-7+(a-l)e 

f 

£ 



In 



0£ + 5(l - 7-x) 
(1 - 7 - a;) [eor(a;) + 5] 



dx ~ -7- In I — 1 + - ^ 



Finally we compute the remaining part of (|92|) 



In 



0£+(l- 7 )5 



(1 - 7)5 + (1 - 7)£<T(i_ 7 )7V +a _i 
0£+ (1-7)5 



1 , 

+ - In 

2 



5 



+ (l-a)5 



5 



+ - In 



2 L(l-7)5 + (l- 7 )^(0)£, 

1 07 -p l(l-a)(7A 2 + 7-4 7 A-A 2 + 2A) 

-£ — — r^; £ = o(l). 



(1 — a)S + (1 — a) ea 
+ 5(1 - a) 



2 |_(1 -«)5+ (1 -a)e<r[l -7 + (a- 1) e] 



2(1-7)75 2 



+ (1 - a) 5] 
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Hence, 



i"[W)]~§M5) + 



2<f>-p-\ 
S 



N -> oo. 



(96) 



Combining (f30"|) . and f3fi]) we obtain 

Vi(S)V2(S)V 3 {S) ~ P(S; iV) 

as in ([85]). ■ 

Next, we shall analyze the function hk($) for N — > oo and (y, z) close to the corner (0, 7). 
To do so, we introduce the new variable I, defined by 



/ = k — c + a, —00 < / < 00, I — 0(1) 
and the new function h^\l, S) = hi +c ^ a (S N) . From |T2|) and (J55|) we have 

(z ' s) = hj exp [r/ K 5iV; 7)] 

c 

Using (J3T|) we see that as iV — > 00 the branch points are located at 

SN - — 

R^SN) ' R 2 {SN) SN 

and that, for iV — > 00 with S and w fixed, 



(97) 



(98) 



exp [r/ (w, SN, 7)] ~ exp 
x exp 



[1- 1 )N + 



S 



In 



SN 



+ 



(l-A)(l- 7 ) 
5 



(99) 



A(l-7) 7 



5w S 

Therefore, we replace the contour C in (J38|) by a Hankel's loop and use (J9T?|) to obtain 



/i (1) (Z,S) 



exp 



1-7 



5 



Ae 



x — </ 
2m J 
n 



exp 



A (1-7) 7 
w 

iSio S 



(l-A)(l-7) 
5 



where 7i begins at w — —00, encircles w = once in a positive sense and returns to its 
starting point. We assume that the branch of U7s" + ^° _1 takes its principal value at the point 
where the contour crosses the positive real axis and is continuous everywhere, except for w 
real and negative 



Introducing the new variable t — , with 



/? = 2VA 7 (l-7) 



(100) 
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we get 



h^(l,S) 



exp 



2S V t 



where 7i' is another Hankel contour of the same shape as 7i. Using an integral representation 
of the Bessel function J.(-) jS2] we obtain 



h k {&) = bP\l,S) 



\+l-a 



(101) 



Using (jHBj) and f)101[) we obtain the following theorem. 



Theorem 5 Let k = c — a + l = \c\ +1, — oo < / < oo, I — O(l) and x — x = 0(1). 
T/ien Ffc(a;) ~ F^ix) as N —>■ oo where 



P 



v^iVY 70(1-7) V27 



l-a 



exp [N$(y)] 



(102) 



Bi^ 



x exp <( — In 



Br + is a vertical contour in the complex plane with Re(S) > 0, J.(-) zs i/ie Bessel function, 
T (•) zs i/ie Gamma function , 

/3 = 2v/A 7 (l-7) 



and 



$(z) = _ z ln(^) - (1 - 2) In (1 - z) + z In (A) - In (A + 1) . 



(103) 



6 The transition layer at y = Yq(z) (Region II) 

As we observed previously in Section HJ the approximation (jUHj) is not valid along the curve 
V = Y o(z) given by 



YJz) 



>> ln( zX + Z ~ y \. 7 <*<1. 



A + l (A + l) 



P 
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The condition y rs Y (z) corresponds to the saddle 9 being close to the pole at i? = in 
jSnj)- To obtain the expansion of (jSHj) for ?/ rs ^(z), we first approximate the integrand in 
f)67)l for small using 



X 2 (l-z) 



7](W. } <&, z) + fi (0) ~ In [(1 - zf- 1 \ z ~ x z~ z ] - Y (z)$ + 



2( 



(A + l) 4 



In 



z \ + z -\\ (z- 7 ) [2pC + 3(A + l)C(^-7) + (A-l)(A+l) 2 (z- 7 ) 2 ] 



(zA + z-A) 2 (A + l)^ 



C = A 2 ( 7 - 1) + 2A - 7 . 



(104) 
(105) 



Using the above in (ffiTj) yields 
1 



F k (x) 



1 1 

exp[$(^)iV] 



1 



2vri y $ 

Br+ 



V2tcN]] z(1-z) 
with given by (JTO3J). Changing variables to 



cxp 



N > dx9 



x) 



V(y,z) = [y-Y (z)} 



N 



Yoz) 



(106) 



we obtain [T2j 



^ | exp|[y-FoW]^iV + ^M#iv|i^ 

Br+ 



V 



2ni 



exp wV + 



v 2 \ 1 



2/w 



exp 



dt 



Br' 



where Br' is another vertical contour, on which Re(v) > 0. Therefore, 

v 



F k (x) = F®(V, 



1 



z) ~ 



y/2nN\ z(l-z 



' exp [$(«)iV] ' 



cxp 



dt. (107) 



This is valid for y — Yq(z) = O (^N 2^ with 7 < 2 < 1 and gives the transition between 
regions !SHi and 9^2 in Theorem El 
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7 The boundary z = (Region III) 



We shall next analyze the part of d® corresponding to z = 0. We start by proving the 
following lemma. 



Lemma 6 Let k = 0(1). Then as N — > oo 

h k (*) ~ ^ (A - 7 ^) fc [-i?i(tf)] y( ^ [i? 2 (0)] [1 - yW]iV , 6+ (0, 0) < i? < 9 

k\ 

where + ($/, z) is defined in 1%5)) . 

Proof. From (|T3j) we have for N — > oo, with w = O(N), 



H($,w) = w N [-^(tf)]^ [# 2 (0)] [1 - yW]Jv exp 



(A - 70) 



iV + O 



iV 

W7 2 



We replace C by C = and change variables from w to w = ^ in (fT2j) to get 



' f IH>>.\<>) 
2m J 



(108) 



(109) 



Using (jl09j) we have 



/i fc (tf)~iV fc [-i?!^]^ [,R 2 



[l-V(i?)]JV 



2vrz 



exp 



(A - 7 0) 



i/ 



i-fe 



■it 



which evaluates to give 
Writing (jHEJ) as 



F fc (x)-F fc (oo) 



A 



N 



X + 1J 2ni 



r) 



Br - 



JV-LcJ-1 



n 



J=0 J 



h k (d)d§ 



and using (}38|) and (|108|l we have 



F k x -F,(oo ~ — 



N 



k\ V A + 1/ 2niJ<dy$-9 

Br" 







(A - 7*)' 



x exp {ytfiV + //(0) JV + V(0) In [-^ (0)] N + [1 - V(#)] In [^(0)] N} d-d. 
Use of the saddle point method yields 



F k (x) - F k (oo) = F^(y) - F fe (oo 



N k ( A 



fc! VA + iy v / 2TiV 



iV 



e V e - e 



(A- 7 e ) fc exp[iVvI/ ( 2/ ,eo)] 



d-d 2 



(y, e c 



;no) 
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where 

(y, ti)=y# + fji{-&) + V(&) In [-RM] + [1 - V(&)] In [R 2 {#)} 
and Qo(y) is the solution to the equation 

d* 



(y, ©o) = o. 



in) 



We note from (jlll|) that 



e (y)^e+(o,o), y ^o 

©o(y) ~> #o, 2/ -> oo 



and thus we shall limit 6 to the range 6 + (0,0) < 6 < 6q in (jllOj) . Expansion (jllUj) 
remains valid for y — >■ and can be used to evaluate Fk(0) — Fk(oo) for k = 0(1). 

8 The boundary z = 1 

We examine the scale k = N — 0(1), which corresponds to 1 — z = 0(iV _1 ). The analysis is 
different for three ranges of y. 

8.1 The boundary layer at z = 1, < ?y < Yq(1) (Region IV) 
Lemma 7 Let k = N - j, j = 0(1), # > 0, & = 0(1). Then, as N -> oo 



l + (l-7)t? 
A 



(112) 



Proof. From (|T2*|) we have 



1 fH($,w) 



2m J wi +1 
c 



dw. 



Changing variables to u — wN we get 

h k (&) = N^Jexp {nV($) In [l - Ri(#)jj] } 
c 



x exp jiV [1 - V(ti)] In l + i2 2 (tf)-^ } 



iV- 7 ® exp 

2m J 

a 



1 + (I- 7R 
A 



-it 



iV -> oo 



where C = i is a small loop around it = 0. Hence, the result follows by expanding 



exp 



l+(l-7)i? 
A 



U 



in Taylor series. 
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Using (HH and (fTT2j) in flU we have 



FJx) ~ — 



JV 



j\ \X + 1J 2niJ$y<d-9 

Br+ 



-00 



i + (i 

A 



exp[y#N + fi(&)N] d-d. 



We evaluate the integral above for N — > oo by the saddle point method and obtain 

'1 + (1-7)6!^ 



F k {x) = Ff\y) 



N> 1 



A \ N 1 



j\ V5kN \x + ij ©i V ©i-^o 

x exp [ye x JV + /i(e x )JV] ' 



A 



where ©i(y) is the solution to 



y + //(eo = o. 



(113) 



(114) 



8.2 The corner layer at (Yo(l), 1) (Region V) 

From ()114j) we see that Qi(y) — > as y — ► lo(l). Therefore, we shall find a solution valid 
in the neighborhood of the point (io(l), 1) . From ()112|) we have as N oo, with j = 0(1) 



1 fN\ j 



and from we get 



7? + 



p[ln(p)-l] + l 

(A + l) 2 L 2(A + lf 

-Y (l)tf+~y 2 (l)tf 2 , ^^0. 



(1- 7 )(P-4A + 1) Cln(p) 



(A + If 



(115) 



A 



1 



Using (J38j) and the above in (J35|) we obtain 

F k (x 



Defining 



A + 1 / 2m J ^ 

Br+ 



exp 



y$N - Y (1)W + W 2 (l)d 2 N 



d-d. 



V(y, l) = [y- Y (l)\ 



N 



Y2O) 



we obtain, for j = 0(1) and y — Yq(1) = O (^N 2 



F k (x)=Ff(V) 



1 / N 



A 



jl \ X J \X + lj v^F 



V(y,l) 



exp 



(it. 



(116) 
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8.3 The boundary layer at z = 1, y > Yq(1) (Region VI) 

We shall examine the portion of the boundary z — 1 for y > Y (l). To do so, we write Fk(x) 
as 

'jv-Lcj-i 



FJx) -F k (oo) 



A \ N 1 



A + 1 J 2m 



Br" 



n 



j=0 3 



h k {§)d$ 



and, repeating the same calculation done in Section IHTTl we conclude that 



F,(x) 



X 



+ A 



N 



N 



(117) 



9 The boundary y = 

In this Section we analyze the part of the boundary where the boundary condition © applies. 
Theorem El is not valid for small y since for 7 < z < 1 the saddle Q(y, z) — > 00 as y — > + , 
but our approximation to the integrand in (JHfij) was only valid for 1? = 0(1). We will scale 
y = x/N and use Lemma EJ which applies for = SN = 0(N). We note that yd = xS. We 
will need to examine separately the cases 7 < z < 1 and 1 — z = 0(e). 



9.1 The boundary layer at y = 0, 7 < z < 1 (Region VII) 
Lemma 8 Let k = Nz, 7 < z < 1, fi = SN, S > 0, 5 = 0(1). TTten, 



y/2irN 



exp < (1 — 2) In 



(*-7)S 

La(i-*) 



AT + (1 - 7) In I - — 1 I A" 



z — 7 



X 



(1-7) 



(1-z) (2-7) 



exp 



_ h) , *-7^ , (i-«)(i-a) 

S \ z — 7 



5 



'118) 



Proof. Changing variables in (|12J) from w to u = wSN, we have 



h k (ti) = / exp (v^SW) In [l - i^SW) 

Z7T2 J I 
C 



'SiVJ 



N-(l-z)\n 



u 

sn 



N 



x exp { [1 - V(SN)} In 1 + R 2 (SN) 



u 
'SN 



A" > ^ 



2vn 



x exp 



exp 



'1 - ^) In 



■it 



iV+(l- 7 )ln(l + ^) JV 



5 V A/ (A + m)5 



U 
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Use of the saddle point method gives, for N — > oo with 5* and z fixed, 



A 



V2nN ui 
x exp 



exp 



(1 - z) In N + (1 - 7 ) In (l + ^) N 



0, (* , «i\ , »! (1 - 7 ) (1 - A) 
^ ln v 1 + TJ + (A + Ul )5 



( l- 7 ) 

(z- 7 ) 3 



where the saddle point is located at 



A(l-^) 

z — 7 



Using (jEEl) and (ITTHj) in (|3E1). we get 



Fu(x) 



1 



A 



2ttA^ V A + 1 



N 



<Pl (1 7) 



exp [—7 ln(7)AT] 



exp jxSAT - (1 - 7) ln[(z-7) SN]N + (1- z)\n 



Br' 



x exp { — In 



7 



_0(z-7)AT 



+ 



[X(l-z) 
p-l + (l-z)(l-A) 



N 



N 



S 



s 



dS. 



(119) 



where Br' is a vertical contour on which Re(S') > 0. Computing the integral in (|119|) as 
N — > 00 by the saddle point method, we find that the integrand has a saddle at 



S*(x, z) 



z — 7 

x 



(120) 



and hence the leading term for ()119j) is 

F k ( X ) r 



1 \ ^ / A ^ 



2ttA^ / V A + 1 



07 (1 — 2) 2 — 7 



exp [—7 111(7) iV] 



x exp < (2 — 7) N — (1 — 7) In 



N 



x 



N + (1 - z) In 



(^~7) S 
A (1 — 0) x 



-iV 



A' 



£0 

x exp < In 

z — 7 



7 



_(p(z-j)N 



20-p-l + (l-^(l-A) 

H .T 



xT ( — — + 1 -a 

2 — 7 



z — 7 

X0 \ 



2 — 7 
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After some simplification we get 



X(j) 

(1 — z) z — 7 \z — 7 



+ 1 - a 



X(j) 

z — 7 



x exp < [z — 7) In 



x exp 



{z-^fN 

X(j) 



N + z\n(\)N - (1 - z) In (1 - z) N - ln(A + 1)JV - 7111(7) JV 

(121) 



In 



2 — 7 



7 



(z - 7) iV 



+ 



2A(l- 7 ) 
z — 7 



+ (a-i; 



Note that (j!21j) is singular as z — > 7. We can show that if we expand (j!21j) as (x, z^j * 
(0,7), with fixed, the result asymptotically matches to he corner layer approximation in 

Theorem |S| The latter must be expanded for x, / — > 00 with - fixed. 
9.2 The corner layer at (0, 1) (Region VIII) 

When z — ► 1, (j!21|) becomes singular. Therefore, we shall find an approximation to Fk(x) 
close to the point (0, 1). We first observe that when 

k = N-j } $ = SN, S>0, j, S = 0(1) 

we obtain from ()12j) . by a calculation similar to that in the proof of Lemma [3 



h k (ti) ~ N 2j - 



(l~7)g 
A 



(122) 



Using 455|) and (IT22I) in Qgft we have 



A' 



2tu 



Br 4 



S \S 



y/2irN\ 07(1-7) 
exp {xSN - (1 - 7) In [(1 - 7) SN] N - 7 ln( 7 )iV} 

;i-7)^ 



(123) 



x exp <( — In 



7 



A 



Using the saddle point method we find that the integrand of (|123|) has a saddle at 

1-7 



S*(x,l) 



X 
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where S*(x,z) was defined in ()120|) . Hence, to leading order, we obtain 



x exp |(1 — 7) In 




P 1 



ex 



1-7 



r 



L(l-7) N\ 

X(f) 



1-7 



+ 1 - a 



2nN V <^y 1 - 7 
iV-7ln(7)AT 

(1-7) S 



(124) 



20 

x exp < In 



7 



Ax 

+ (3A-l)x}>. 



1-7 [0(1-7)A^_ 
Note that from ((EH) and (fl24jl we have 

F fc (x) = O (x Q+(z ~ 7)7V ) =0(x fc - [cJ ), x^O, LcJ + 1 < A; < 

and thus we recover the result found in [2j (equation (38) therein), that 

d n F k 



dx n 



.(0) = 0, [c\ + 1 +n < k < N. 



10 Summary and numerical studies 

In most of the strip D = {(y, z) : y > 0, < 2: < 1} , the asymptotic expansion of F^ix) is 
given by Theorem El Below we summarize our results for the various boundary, corner and 
transition layer regions, where Theorem El does not apply. The paragraph number refers to 
the corresponding region (see Figure EJ. 

Theorem 9 As N — > 00, asymptotic approximations of the joint distribution F^{x) are as 
follows. 

I. k — c — a + I — [c\ + I, —00 < / < 00, I — 0(1) and x — j^, X — 0(1), 
Fk(x) ~ Fi(x) as N — > 00 where 



^ (1) (x) 



l-a 



5 V 5 



Br 4 



x exp <^ — In 







/3S 



2A0(1- T ) 



5 I V5 



(125) 
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0.2 



0.4 



0.6 



0.8 



Figure 2: A sketch of the different asymptotic regions. 

Br + is a vertical contour in the complex plane with Re(S') > 0, J.(-) is the Bessel 
function, Y (•) is the Gamma function , 

/3 = 2VA 7 (l-7) 

and 

$0) = -z ln(z) - (1 - z) In (1 - z) + z In (A) - In (A + f ) . (f 26) 

//. y - Y Q (z) = O (n~*) , 7 < z < 1 



«'» = ~ ^V^j exp| * WjV| ^ / exp (-0 * (127) 



-oo 
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yJ A + l (A + l) 2 V P 



, 7 < z < 1 



Y 2 {z) 



2C 



2 — 7 



(a + i)(\z + z - xy 



(A + l) 
2CP 



In 



z + zX — X 
P 



+ J^(z- 1 ) + (X-l)(z- 1 ) 2 



(A + l) 2 ' (A + l) 
C = 2A- 7 + (7-l)A 2 . 



///. k = 0(1) 



N k / \ \ N 1 



X- 



1 



-e 



e V e - e ( 



(A- 7 e ) fe exp [jV* (y, ©o)] 



^2 



(1/, ©0) 



where 



^0 (l/, = + + In (0)] + [1 - V(&)] In [^(0)] 
and 60 (y) is the solution to the equation 



(y, ©o) = 0. 



IV. k = N-j, j = 0(1), 0<y<Y (l) 



F k (x) = Ff(y) 



Ni 1 



X \ N 1 



-0n 



j\ V2^N \X + 1J 0i V ©i -^0 
x exp [yQiN + //(6i)iV] ' 



1 + (1-7)0! 
A 



VpWi) 



where Q\(y) is the solution to 



y + ^'(Q 1 )=0. 
V.k = N-j, 3=0(1), y - Y (l) = O (n~^ 



F k (x) = Ff\v) 



1 (N 



X 



j\ v XJ \X + lJ 



JV 



V(y,l) 



ex P I — g 



with 



V(y,l) = [y-Y (l)} 



N 
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VI. k = N-j, j =0(1), y>Y (l) 



F k (x)-F h {oo) = Ff>(y) 



1 + A 



iV 



N 



(133) 



where F^ (y) is as in item IV. 
VII. x = 0(1), 7 < ^ < 1 

F fc (x) =F^(x,z) 

x exp < (z — 7) In 



27riV 



r 



07 (1 — z) z — 7 V- 2- 7 



£0 



+ 1 - a 



X(j) 

z — 7 



x exp 



_(;z-7) 2 iV 

£0 



N + z\n(\)N - (1 - z) In (1 - z) N - ln(A + l)N - 7ln( 7 )iV 

(134) 



In 



2 — 7 



7 



(2 - 7) N 



2A(l- 7 ) 

2; — 7 



+ (A - 1) 



x 



VIII. k = N-j, j = 0(1), x = 0(1) 



1 / A 



N 



x exp < (1 — 7) In 



P 1 

j! VA + I/ 2vriV\/ 07 1 - 7 
iV- 7 ln(7)iV 



(1-7) N. 



(135) 



x(f) \ a f X(f) 



1-7 



x exp 



£0 



r 



In 



1-7 



+ 1 - a 



7 



g-7r 

Ax 

+ (3A-l)x}>. 



1-7 [0(1-7)A^_ 
From (j3T?j) we see that the continuous part of the density (for x > 0) is given by 



f k {x) = F' k {x) 



A \ N 1 



A + 1 7 2m 



Br 4 



II 



0, - i? 



h k (d)dd. 



(136) 



Therefore, we can approximate the density by multiplying the asymptotic approximations 
by the corresponding saddle point: 

f k (x)~QG x (y,z) 
f k {x)~e+G 2 (y,z) 

f k (x)~O [Fi 3 \y)-F k (oo) 

( 4 )/ 



f k (x)~e lF ^(y) 
f k (x)~S*(x,z)NFM(x,z). 
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x 



Figure 3: A comparison of the exact (solid curve) and asymptotic (+++) values of the 
density fk{x) for k = 17 with N = 20. 



Here the right sides are the expansions valid in the various regions. 

We compare our asymptotic approximations to the density with the exact results obtained 
by evaluating ffHfij). We use N = 20, 7 = 0.37987897 and A = 0.0122448. In Figure H we 
graph fk(x) and the asymptotic approximation QGi(y,z), for k = 17. We note that the 
asymptotics predict that the maximum value will be achieved at x = NY (J|J ~ 3.052. 

In Figure |U we graph fs(x) and the asymptotic approximation B+G^y, z). We note that 
the maximum is now achieved at x — 0, which is consistent with the asymptotic analysis. 

In Figure El we graph fo(x) and 9 Fq 3 \i/) — F (oo) . This shows the importance of 
treating Region III. 

In Figure El we graph In [/2o(^)] an d In QiFq A \u) , the latter being the boundary layer 
approximation to the density at z — 1. 

In Figure [7|we graph In [/^(x)] and In [S*(x, z)NF^ 7 \x^ z)\ with x = 0.001. We see how 
the asymptotic approximation breaks down for 2^7 and for z ~ 1 (it does not apply for 

z < 7). 
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Figure 5: A comparison of the exact (solid curve) and asymptotic (ooo) values of the density 
f k (x) for k = with N = 20. 
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Figure 6: A comparison of the exact (solid curve) and asymptotic (o o o) values of In [/&(#)] 
for z = 1 with N = 20. 
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Figure 7: A comparison of the exact (o o o) and asymptotic (solid curve) values of In [f k ( 
for x = 0.001 with N = 20. 
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Figure 8: A comparison of the exact (+++) and asymptotic (solid curve) values of the 
density fk(x) for x = 5 with N = 20. 

In|S]we graph fk(x) and the asymptotic approximations QGi(y, z) and Q + G2(y, z) with 
x = 5. To show in more detail the different ranges, in Table 1 we compare the exact value of 
the density with the approximations given by QqF^ (y) and Q + G2{y, z), for x — 1 (y = 0.05) 
with < k < 9. In Table 2 we compare fk(x) with the approximations given by QiF^ 4 \y) and 
QGi(y, z), for x = 1 (y = 0.05) with 10 < k < N. We note that for our choice of parameters 
c ~ 7.598 and Y (l) ~ .4998. Also, when y = 1/N = 0.05, the equation Y*(z) = y has the 
solution z ~ .4811, so that the transition between d\ 2 and occurs at Nz ~ 9.622. 

11 Conditional limit laws 

In this section we will study the probability of k sources being on, given that the buffer 
content is equal to x > and the probability of the buffer having some non-zero content x, 
given that k sources are on, for various ranges of x and k. 
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Table 1: 



k 





Exact 




® + G 2 (y, 


*) 





-3.418 


.445410" 


-19 


.432410" 


-19 






1 


-3.593 


.125510" 


-17 


.113310" 


-17 


.133410" 


17 


2 


-3.785 


.168610" 


•16 


.148610" 


-16 


.172710" 


16 


3 


-3.994 


.142510" 


■15 


.129810" 


-15 


.141510" 


15 


4 


-4.216 


.835510" 


-15 


.850310" 


-15 


.818710" 


15 


5 


-4.444 


.352210" 


-14 


.445910" 


-14 


.350110" 


14 


6 


-4.660 


.106210" 


-13 


.194710" 


-13 


.101010" 


13 


7 


-4.825 


.219610" 


-13 






.207310" 


13 


8 


-4.838 


.283310" 


-13 






.270710" 


13 


9 


-4.411 


.182610" 


-13 






.164810" 


13 



Table 2: 



k 





Exact 


eGi(i/, 


*) 


0i^ 4j (y) 


10 


-3.177 


.238910- 


-14 


.229710- 


-14 




11 


-1.650 


.726610" 


-16 


.740010' 


-16 




12 


-.2175 


.857810" 


-18 


.863510" 


-18 




13 


1.115 


.497210" 


-20 


.507410- 


-20 




14 


2.378 


.159910" 


-22 


.164610" 


-22 




15 


3.591 


.304410" 


-25 


.314510- 


-25 




16 


4.769 


.354010" 


-28 


.374010- 


-28 


.146710- 27 


17 


5.921 


.251910" 


-31 


.265010- 


-31 


.532610" 31 


18 


7.052 


.106310" 


-34 


.114310" 


-34 


.144910" 34 


19 


8.167 


.242210" 


-38 


.242210- 


-38 


.263010~ 38 


20 


9.269 


.228410" 


-42 






.238610" 42 
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11.1 The distribution of active sources conditioned on x 

We define 

M 1 (x) = y j f k (x) = Km — Pi[X(t)<x], x>0, (137) 
^ — ' t^oo ax 

k=0 



M 2 (k) = [ f k {x)dx = lim Pr [Z(t) = k, X(t) > 0] , < k < N (138) 
J 



o 



and the conditional probabilities 

A^ =lim 

Mi (a;) t^oo 



f(k | x) = = lim Pr [Z(t] = k \ X(t) = x > 0] , x > 0, 0<k<N (139) 



fix | k) = = hm ^- Pr [X(t) < x | Z(t) = k, X(t) > 0] , x > 0, < fc < iV. 

(140) 

Note that 

M 2 (k) = F fc (oo) - F fc (0), < fc < iV. (141) 

We first consider a; = O(N). In this range we compute Mi(x) using the approximation 
in 9^3 



£ /fc(;r) ~ 2tt ( A + l) / W + (Q+, z) V 0+ - 9 \] Vww (W + , 0+, z) 







X 



* w (y,W + ,0+,z) 



exp [A^(y, W+,0+,,2)] d*. 



The main contribution to the integral comes from near z = 7 and by Laplace's method we 
have 



1 / a \ N 1 -e a 



X 



1 expjiV^ [y,l,0+(y,7), 7 ]} 



since 



[s/,l,e+(y,7), 7 ] 

^(y,W + ,0+ z) = for z = 1 

d 2 ^ + 
■^- T (y,W + ,0 + ,z) = — for z = 7 
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and 

Therefore, 
f(k | x) 

X . 



W + (#, 7 ) = l, WM,7) 



'e+ (y,i)-e 



s/2^Nw + (e+,z)V e+-e u Vww (w + ,e+,z) 



'^[y,l,e+(y,7), 7 ] 



and expanding this for z — > 7 yields the Gaussian form 

1 /-© + (l/,7) 



exp {iV*(2/, W + , e+, z) - iVv]/ [y, 1, 9+ (y, 7) , 7] } 



/(Jfe I x ) = f [Nj + VNA I Ny 



V2nN 



P 



exp < - 



2 dz 2 



[l/,l,e + (y, 7 ),7] (z--y) 2 N 



1 /-e+(2/,7) 



y/2-rcN 



P 



1 e+(y, 7 ) A 2 



6XP ^ 2 „ 



with A = v AT (2; — 7) . This means that given a buffer size x = Ny, y > 0, the most likely 
number of active sources will be c = A/7, with a Gaussian spread whose variance is obtained 
by solving (JH3J) with z — 7. This is also illustrated in Figure |H1 

Next we assume a small buffer size, with x = 0(1/N). Again the main contribution to the 
sum ()137|) comes from k m c, but for x = 0(1/N) we must use the Region I approximation 
f k (x) ~ NF^'ix)- Thus, 



Z-a 



X 



5FV7*(i-7)V5f>l eXP|jV * (7)l 



2m 



s 



Br 4 



x exp <( — In 



which we can write as a residue series 



ps 



fk(x) 



2A0(l- 7 ) 



l-a 



5 



5 



P< ^ ' tt" ) exp [JV$(t)] 



2tt V 7(1 -7) V 2 7 



x y^exp 

3=0 



j + l-a 



X 



+ a) (p 



(142) 



x 



+ a) 



3-1 



" 2A(l- 7 ) 



j+l-a 



Jl-j-1 



(a-1- j] 



P 
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where we have used 
Re s 

Using the generating function (52] 



r(^ + i-Q , s = — i — ; 

b } a — 1 — j 



(-1) J+1 
j\ (j ■ 1 a) 



• , j>0. 



j{x)z J = exp 



3=-oo 



X ( 1 

— I z 

2 V z 



in (11421) we obtain 



and thus 



Z=-oo V ' 



M 1 (x) 



l-(J+l) 



Jl-j-1 



(a — 1 — j) — 



exp 



iV / p0 



27r V 7 (1 - 7) 

j-i 



(j + 1 - a 
x — exp 



exp [N$(j)] y^exp 

j=0 

2p 



J- 



-(j + l-a) 



j + 1 - a X _ 

A(l- 7 ) 13 " +1 " a 

7 



It follows that 



ft) E ^P 



f(k\x) = f(c + l-a\£) 



± v 1 (j+i-q)(p-<? 



j+l-a 



j! 



j'-i 



2A(l- 7 ) 

/3 



j+l-a 



j-1 



(a-l-j) 



Eexp -t^X+ + 



j=o 

which is a discrete distribution in I. 



(j+i-oQ 



A(l~7) 

7 



j+l-a 



11.2 The buffer density conditioned on k 

We first consider 7 < 2 < 1. In this range i^(0) = and from (|141|) we have M 2 (k) = Fk(oo). 
We have 

^[ y ,W.,Q,z] = Q(y,z) 

so that the maximum of the density occurs at = 0, which corresponds to the transition 
curve 

z — 7 



y = Y Q (z) 



A + l (A + l) 



p 1 zX + z - A 

In ( ) , 7 < z < 1. 
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Therefore, we shall use the approximate solution F^ 2 '(V,z), since Region II corresponds to 
6^0. From (jl27|) we get, for 7 < z < 1 



fk{x) 



1 dF® dV 



N dV dy 2-nNV z(l-z] 



exp [$(z)N] 



I ( V 2 

exp 



Y,(z 



Yo(z) 



2( fz + zX-X 
In 



z — 7 



(A + l) {Xz + z-Xy 
From (j!41|) we obtain 



(A + l) 4 
2Cp 



(A + l) 2 (A + l 



P 

' k ( z -j) + (X- 1) (z - 7 ) 2 



and we conclude that 

/(a: \k) = f NY (z) + ViVfi I Nz 



2ir ^Y 2 {z)N 



exp 



1 Q 2 



2Y 2 (z) 



with f2 = [y — Y (z)] v iV. This shows that if the number of active sources exceeds the service 
rate, the buffer content is likely to be large, approximately NY (z) and with a Gaussian 
spread proportional to Y 2 (z) (see also Figure EJ). While this calculation assumed that 7 < 
z < 1, we can show that the result remains valid for z — > 1. 

Next we consider the range < z < 7. We use the asymptotic approximation Fk(x) — 
F fc (oo) ~ G 2 (y,z), or 



/fc(z) 



A 



N 



1 



-00 



1 



2ttN \x + 1; w+(e+, ^) V e+ - 0o V W^V, © + , *) 
exp [iv*(?/,iy + ,e + , 



From (j!41j) we have 
Hence, 



* M (y,W+,e+,z) 

M 2 (k) = F k (oo) - F fc (0) ~ -G 2 (0, ^ 



/(x 1 fc) ~ -e + (o,z 



w + [Q + (o,z),z] e+(o,z)-e Nf M [o,w+,e+(o,z),z 



w + (e+,z) 



x ^ ; ^ e+(o, z), g] exp {jV ^ ^ ^ @+) z] _ m [Q) @+(0; ^ } 



53 



and as y — > 



f(x | k) = f{x | Nz) ~ -0 + (O,z)exp 
= -0 + (O,z)exp [9 + (0,z)x] . 



iV— [0,W + ,Q + (0,z),z]y 



This shows that if the number of on sources is less than the service rate c and the buffer 
is not empty, the conditional density is approximately exponential in x, with mean _ e+ ^ ^ • 
This exponential shape is also illustrated in Figure 0] 

Finally, we consider the case k pa c. We shall again use the approximation in Region I 
and from ()142|) we get 



M 2 (k) 



l-c 



(j + 1 - a) (p - 



x 



(j + 1 - gO 



2A(l- 7 ) 

/3 



J, 



i-j'-i 



(a - 1 - j; 







Therefore, 



/(x | fc) = / (-^r I c + Z-a 



4>N exp 
i=o 



i- — y 



O'+l-«)(p-0)" 


(j- + l_ a )j-i 


"2A(l- 7 )' 




ii 


/3 



1 i+l-a 



J/ 



E ex P 



"(j+l-a)(p-(/>)" 




"2A(l- 7 )" 


J 




/3 



-1 j+l-Q 



J; 



which is a density in x, consisting of an infinite mixture of exponentials. 
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